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Abstract 


The  Behind- Armour  (BA)  debris  cloud  created  behind  the  target  after  a  projectile 
impact  can  easily  be  observed  using  x-ray  technology.  Some  countries  use  this 
technology  to  evaluate  the  lethality  of  the  threat  since  the  number  of  fragments,  mass 
and  velocity  can  be  estimated.  Usually,  two  orthogonal  x-ray  stations  are  located 
behind  the  target  in  order  to  obtain  orthogonal  images  at  two  different  times.  Since 
only  projected  images  of  BA  debris  are  recorded  on  fdms  (two  dimensions),  a  major 
issue  to  solve  concerns  the  available  information.  An  analysis  of  these  x-ray  films 
becomes  time-consuming  when  the  fragments  must  be  matched  from  one  station  to 
another  or  from  corresponding  horizontal  and  vertical  views.  Also,  the  distances  seen 
on  the  films  are  projected  distances  that  are  different  from  the  real  distance  travelled 
by  fragments.  This  memorandum  presents  equations,  methods  and  tools  that  can  be 
used  to  rapidly  locate  fragments  in  space  (three  dimensions)  in  order  to  accurately 
compute  their  velocity,  trajectory  and  finally  assess  their  mass  based  on  x-ray  films 
only. 


Resume 


Un  nuage  de  debris  genere  derriere  une  cible  a  la  suite  de  Timpact  d’un  projectile  peut 
facilement  etre  observe  en  utilisant  la  technologic  des  radiographies  eclair.  Quelques 
pays  utilisent  cette  technologic  pour  evaluer  la  letalite  d’une  menace,  puisque  le 
nombre  de  fragments,  leur  masse  et  leur  vitesse  peuvent  etre  determines. 
Habituellement,  deux  stations  orthogonales  sont  situees  derriere  la  cible  afm  d’obtenir 
des  images  orthogonales  pour  deux  temps  differents.  Un  probleme  majeur  a  resoudre 
concerne  Tinformation  disponible,  car  seule  Timage  projetee  des  fragments  est 
enregistree  sur  les  films  (deux  dimensions).  Analyser  ces  films  de  radiographic 
necessite  demande  beaucoup  de  temps  lorsqu’on  cherche  a  jumeler  les  fragments 
d’une  station  a  Tautre  ou  du  film  horizontal  au  film  vertical  correspondant.  De  plus, 
les  distances  sur  les  films  sont  celles  projetees  et  non  les  distances  reelles  parcourues 
par  les  fragments.  Ce  memorandum  technique  presente  des  equations,  methodes  et 
outils  qui  peuvent  etre  utilises  pour  localiser  rapidement  les  fragments  dans  Tespace 
(trois  dimensions)  afm  de  calculer  avec  precision  leur  vitesse,  leur  trajectoire  et 
fmalement  determiner  leur  masse  uniquement  a  partir  des  films  de  radiographic  eclair. 
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Executive  summary 


A  series  of  ballistic  tests  were  performed  at  DRDC  Valcartier  (Defence  R&D  Canada 
-  Valcartier)  to  study  Behind-Armour  (BA)  debris.  These  BA  debris  were  recorded  on 
x-ray  fdms  and  witness  packs.  X-ray  films  make  it  possible  to  visualise  the  debris 
cloud  at  a  specific  time.  Using  orthogonal  stations  the  fragments  can  be  observed  from 
two  different  directions.  The  x-ray  heads  emit  rays  that  are  absorbed  by  fragments 
creating  white  areas  of  films.  The  further  the  fragment  is  from  the  film,  the  bigger  the 
projected  fragment  image  will  be.  Also,  the  further  the  fragment  is  toward  the  left  or 
right  of  the  emitting  head,  the  more  will  the  projected  distance  be  different  from  the 
real  distance.  Only  these  projected  area  and  distance  are  available  at  the  end  of  the 
tests. 

To  summarise,  to  assess  the  fragment  velocity  or  trajectory  with  the  use  of  the  x-ray 
films,  one  has  to  mathematically  transform  the  projected  distance  to  find  the  real 
distance  travelled.  This  is  very  tricky  since  the  fragments  travel  in  space  (3 
dimensions)  and  the  user  must  deal  with  two  orthogonal  projection  views  (2 
dimensions)  to  define  the  trajectory.  Therefore,  up  to  now  at  DRDC  Valcartier,  only 
the  projected  distance  on  a  film  was  used  to  calculate  the  fragments  velocity.  Matching 
fragments  for  one  station  to  another  and  from  vertical  to  horizontal  film  was 
performed  without  electronic  devices,  by  looking  only  at  the  fragment  shape  and  the 
trajectory  on  films.  Obviously,  this  was  time-consuming  and  errors  could  easily  be 
made. 

This  memorandum  present  equations,  methods  or  tools  to  match  fragments  from  a  film 
to  another,  to  evaluate  the  fragment’s  trajectory  and  velocity  (using  the  real  distance 
travelled),  and  finally  to  estimate  the  fragment’s  mass.  This  can  be  done  since  the 
fragment’s  position  in  space  will  be  well  defined. 

Since  witness  packs  are  often  used  in  ballistic  tests  to  record  BA  debris  and  assess 
fragments  mass  and  velocity,  the  information  obtained  with  x-ray  can  now  be  used  to 
check  or  validate  the  results. 

Normal  impact  tests  using  x-ray  were  performed  at  DRDC  Valcartier  to  validate 
equations  and  models  presented  in  this  memorandum.  Moreover,  an  Excel  program 
was  built  in  order  to  reduce  the  analysis  time. 


Y.  Baillargeon  and  C.  Lalanne.  (2005).  Methods  to  perform  behind-armour  debris 
analysis  with  x-ray  films.  DRDC  Valcartier  TM-2003-123.  Defence  R&D 
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Sommaire 


Plusieurs  essais  balistiques  ont  ete  effectues  a  RDDC  Valcartier  (R&D  pour  la  Defense 
Canada  -Valcartier)  pour  etudier  les  debris  generes  derriere  la  cible.  Ces  debris  ont  etc 
enregistres  sur  des  fdms  de  radiographic  eclair  et  sur  des  panneaux  temoin.  Les  films 
de  radiographic  permettent  une  visualisation  du  nuage  de  debris  a  un  temps  precis.  En 
prenant  des  stations  orthogonales  les  fragments  peuvent  etre  observes  a  partir  de  deux 
directions  differentes.  Quand  les  fragments  passent  devant  les  films,  les  tetes  de 
radiographic  emettent  des  rayons  qui  sont  absorbes  par  ceux-ci  creant  ainsi  une  zone 
blanche  sur  le  film.  Plus  le  fragment  est  loin  du  film  plus  I’image  projetee  de  celui-ci 
sera  grosse.  Par  ailleurs,  plus  le  fragment  est  a  gauche  ou  a  droite  de  la  tete  emettrice 
et  loin  du  film,  plus  les  distances  projetees  lues  sur  le  film  seront  loin  des  distances 
reelles.  Seules  les  informations  sur  les  aires  et  les  distances  projetees  sont  disponibles 
a  la  fin  des  tests. 

En  resume,  pour  calculer  la  vitesse  ou  la  trajectoire  des  fragments  en  utilisant  les  films 
de  radiographic  eclair,  on  doit  transformer  mathematiquement  les  distances  projetees 
pour  trouver  les  distances  reellement  parcourues.  Ceci  est  vraiment  difficile  a  faire 
puisque  les  fragments  voyagent  dans  I’espace  (3  dimensions)  et  I’utilisateur  doit 
travailler  avec  deux  vues  qui  sont  des  projections  orthogonales  afin  de  reconstituer  le 
trajet.  Par  consequent,  jusqu’a  maintenant  a  RDDC  Valcartier,  le  calcul  de  la  vitesse 
des  fragments  etait  effectue  en  utilisant  seulement  la  distance  projetee  sur  un  des  films. 
Ee  jumelage  des  fragments  d’une  station  a  une  autre  et  du  film  horizontal  au  film 
vertical  etait  fait  sans  I’aide  d’outils  electroniques,  c’est-a-dire  uniquement  en 
regardant  la  forme  du  fragment  et  la  trajectoire  parcourue  sur  les  films.  Evidemment 
ceci  necessitait  beaucoup  de  temps  et  des  erreurs  pouvaient  facilement  se  produire. 

Ce  memorandum  presente  des  equations,  methodes  ou  outils  pour  permettre  le 
jumelage  de  fragments  d’un  film  a  un  autre,  pour  evaluer  la  trajectoire  et  la  vitesse  de 
fragment  (en  utilisant  la  distance  reelle  parcourue)  et  finalement,  de  determiner  la 
masse  des  fragments.  Ce  sera  possible,  car  la  position  du  fragment  dans  I’espace  sera 
bien  definie. 

Puisque  les  panneaux  temoins  sont  souvent  utilises  dans  les  essais  balistiques  pour 
enregistrer  les  debris  generes  derriere  la  cible  et  determiner  la  masse  et  la  vitesse  de 
ceux-ci,  I’information  obtenue  a  partir  des  radiographies  eclair  peut  maintenant  etre 
utilisee  pour  verifier  ou  valider  les  resultats. 

Des  tests  a  angle  d’impact  normal  ont  etc  effectues  a  RDDC  Valcartier  a  I’aide  de 
radiographic  eclair  afin  de  valider  les  equations  et  modeles  presentes  dans  ce 
memorandum.  De  plus,  un  programme  Excel  a  etc  cree  pour  reduire  le  temps 
d’ analyse. 


IV 
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1.  Introduction 


When  a  projectile  impacts  a  target  in  experimental  trials,  a  system  is  required  to  record 
the  ballistic  characteristics  of  this  event  before  impact  and  behind  the  target.  Three  of 
those  characteristics  are  the  projectile  velocity,  integrity  and  yaw.  X-ray  technology  is 
often  used  to  accurately  obtain  all  this  information.  Indeed,  x-ray  photographs  taken  at 
different  times  give  information  on  the  projectile  flight.  DRDC  Valcartier’s  Small 
Caliber  Laboratory  has  a  system  that  is  composed  of  two  orthogonal  x-ray  stations  and 
fdms  located  before  and  after  the  target.  Therefore,  using  this  system  it  is  possible  to 
find  the  real  projectile  velocity  in  a  semi-automatic  way  knowing  the  distance  between 
the  x-ray  stations  and  the  time  required  by  the  projectile  to  go  from  the  first  station  to 
the  other.  The  total  yaw  can  be  obtained  with  the  orthogonal  views  and  the  integrity 
check  on  all  films. 

However,  some  issues  must  be  solved  in  order  to  assess  individual  Behind-Armour 
debris  velocity  easily.  This  memorandum  presents  mathematical  equations  to  solve 
those  issues.  The  first  issue  is  to  correctly  associate  fragments  on  horizontal  and 
vertical  films.  Next,  the  association  of  each  fragment  from  one  station  to  the  other 
must  be  made.  The  impact  point  on  the  witness  plate  also  needs  to  be  found  in  order  to 
associate  holes  in  a  witness  pack  to  the  x-ray  films.  Once  these  issues  have  been 
solved  and  the  fragment  trajectory  defined,  the  next  step  is  to  compute  real  fragment 
velocity. 

Of  equal  concern  is  to  estimate  fragment  mass  with  the  available  information.  The 
photographs  give  the  orthogonal  fragment  views  at  different  times.  Therefore,  four 
images  are  available  with  the  hole  area  created  by  the  projectile  on  the  witness  plate  to 
assess  the  mass.  Equations  are  developed  in  this  memorandum  that  will  allow 
fragment  characterisation  (mass,  velocity,  trajectory). 

Software  have  already  been  developed  that  use  similar  approach  to  analyse  x-ray  [1,2]. 
However,  other  techniques  exist  to  obtain  debris  cloud  characteristics  and  all  present 
advantages  and  constraints.  Amongst  these  are  the  holographic  approach  [3],  switch 
screen  approach  [4]  and  finally  an  approach  using  a  laser  beam  and  CCD  cameras  [5]. 
The  last  approach  uses  a  similar  technique  to  the  x-ray  method  described  here. 
However,  this  memorandum  will  be  limited  to  the  x-ray  issues  to  solve. 
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2.  Set-Up 


Diagrams  of  the  x-ray  configuration  found  in  the  DRDC  Valcartier  Small  Caliber 
Laboratory  are  presented  in  Figure  1,  2  and  3.  The  projectile  comes  from  the  left  and 
impacts  the  target.  This  results  in  the  ejections  of  multiple  fragments  behind  the  target. 
When  these  fragments  are  in  front  of  the  films,  x-ray  heads  of  orthogonal  stations  1 
and  2  are  activated.  Finally,  fragments  end  their  trajectory  on  the  witness  plate  (WP) 
on  the  far  right.  The  parameters  indicated  on  Figure  1,2  and  3  are  presented  in  Table  1. 
These  parameters  are  constant. 
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Figure  1.  Side  view  at  normai  impact 
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Figure  2.  Upper  view  at  normai  impact 
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Table  1.  Set-up  parameters 


PARAMETER 

SIGNIFICATION 

a 

Vertical  distance  between  the  x-ray  station  and  the  line  of  fire 

g 

Vertical  distance  between  the  line  of  fire  and  the  horizontal  film 

C 

Vertical  distance  between  the  line  of  fire  and  the  vertical  film 

e 

Horizontal  distance  between  the  horizontal  x-ray  station  and  the  line  of  fire 

b 

Horizontal  distance  between  the  impact  point  on  target  and  the  reference  origin  (x 
distance) 

i 

Horizontal  distance  between  the  two  X-ray  stations 

j 

Horizontal  distance  between  the  impact  point  behind  the  target  and  the  WP  center 

9  (theta) 

Target  degree  of  inclination 

Figures  4  and  5  present  the  horizontal  and  the  vertieal  radiography  pietures  of  behind- 
armour  debris.  Eaeh  figure  presents  the  double  exposure  of  the  residual  penetrator  and 
behind-armour  debris  at  two  different  times.  An  undeformed  residual  penetrator  and 
pieces  of  the  jacket  impacted  the  WP.  In  this  case,  the  jacket  is  completely  stripped 
from  the  core.  The  forward  segment  of  the  jacket  is  broken-up,  whereas  the  tail 
segment  is  almost  intact. 
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Figure  4.  Horizontal  x-ray  view  0.5  APM2  ammunition 
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Figure  5.  Vertical  x-ray  view  0.5  APM2  ammunition 
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3.  Fragment  projections  on  films 


The  first  issue  to  be  solved  is  the  association  of  the  fragment  projection  on  films. 
Figure  6  presents  a  3D  representation  of  the  location  (PI,  P2,  P3)  of  a  fragment  at 
different  times  relative  to  the  x-ray  heads.  Indeed,  the  fragment  PI  is  projected  on  two 
films,  one  placed  horizontally  (H)  and  the  other  vertically  (V). 


The  characteristic  that  distinguishes  the  fragment  projected  images  Ply  and  Pin  from 

any  other  fragment  images  on  the  film  is  that  their  vectors  ( D\  and  D2 )  are  crossing 
at  point  PI.  These  two  vectors  are  made  of  a  line  joining  the  fragment  image  on  the 
film  and  the  x-ray  source.  Equation  1  is  used  to  evaluate  the  smallest  distance  between 

vectors  D\  and  D2  and  will  be  used  as  a  criterion  to  confirm  that  the  corresponding 
fragment  projected  images  are  matched  together. 

A.,i.  =  pK-‘i2,f  +(rfi, -d2,:f  +(dL -d2:f  (1) 

In  general,  the  smaller  the  distance,  the  better  is  the  match.  The  fragment  positions  on 
the  films  were  given  according  to  the  reference  coordinate  system  (x,  y,  z)  located  at 
the  virtual  interception  point  of  the  line  of  fire  and  the  normal  to  both  horizontal  and 
vertical  X-ray  heads  front  lens  at  the  first  orthogonal  x-ray  station  (see  Figures  1,  2  and 

3).  Equations  2  and  3  present  the  equations  of  vectors  D\  and  D2  that  need  to  be 
solved  to  define  the  coordinates  of  the  point  on  each  vector  that  is  the  closest  to  the 
other  vector.  The  coordinate  of  this  point  might  not  be  exactly  the  same  for  both 
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vectors  (i.e.  leading  to  a  value  of  Dmin  equal  to  0)  even  if  the  right  match  has  been  done 
since  there  is  a  tolerance  in  measurement. 


D\  =  R,,+t^V\ 


dl, 

d\^ 


(2) 


D2  =  R,, 


+  sV2 


d2^ 

d2^ 

d2^ 


(3) 


Where  D\  and  D2  give  the  parametric  equations  of  a  line,  V\  and  V2  are  the 
unitary  vectors  and  t^  and  Sm  are  the  scalar  value  that  allow  positioning  on  the  lines. 
Ry^  and  R^^  are  presented  in  Equation  4  and  5  and  represent  the  x-ray  stations’ 
coordinate  in  space. 


R 


VI 


0 

0 

a 


(4) 


Where  Ry^  is  the  first  vertical  x-ray  station  position  on  the  side  view  (Figure  2)  and 
a  is  the  distance  between  the  x-ray  station  and  the  fire  line. 


R 


H\ 


(5) 


Where  is  the  first  horizontal  x-ray  station  position  on  the  upper  view  (see 
figures)  and  e  is  the  distance  between  the  horizontal  x-ray  station  and  the  fire  line. 

Equations  6  and  7  are  used  to  compute  the  unity’s  vectors,  VI  and  V2  . 


Xfh 

y  FH 

_-(g  +  a)\ 

'J^FH  y  FH 

(6) 
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V2 


Where  V\  and  V2  are  the  unitary  vector  required  in  the  parametric  Equations  2  and 
3,  respectively,  and  Xpy,Zpy  are  parameters  giving  the  fragment  centroid 

position  on  the  horizontal  and  vertical  film,  respectively,  and  finally,  parameters  g,  a, 
c  and  e  were  defined  in  table  1 . 

In  order  to  solve  Equations  2  and  3,  the  two  last  parameters  and  must  be  found. 

The  distance  between  D\  and  D2  is  minimum  when  the  gradient  of  D\  —  D2  is 
equals  to  zero.  Equations  8  and  9  present  the  solution  in  this  case. 

J  ^  ^  -R^,)V\  +  {VUV2){Ry,  -R^,)V2\  (8) 

\-{V\»V2f 

Where  is  a  scalar  value  for  the  parametric  Equation  2. 

J  \Ry,-R,,)T2-{VuT2){Ry,-R,,)V^  (9) 

l-(Vl*V2y 

Where  is  a  scalar  value  for  the  parametric  Equation  3. 

Therefore,  it  is  possible  to  associate  the  horizontal  fragment  projected  image  that  best 
match  with  the  vertical  fragment  projected  image  using  all  these  equations. 

Once  the  best  match  has  been  defined,  the  real  fragment  position  in  space  must  be 
estimated  using  Equation  1 0  since,  as  mentioned  previously,  equations  2  and  3  will 
not  give  the  exact  same  coordinate.  To  evaluate  the  best  real  fragment  position  (PI), 
the  middle  of  the  virtual  line  having  [dfr,  dly,  dfr]  and  [d2x,  d2y,  d2z]  as  extremities 
is  selected. 

Ud\,+d2^)l2' 

P\=  =  {d\y+d2y)l2  (10) 

zj  [(r/l^+J2J/2 
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4.  Fragment  association  between  stations 


Another  issue  is  to  assoeiate  a  fragment  from  an  x-ray  station  to  another.  As  a  first 
step,  projected  fragment  images  should  have  been  matched  by  pairs  for  each  station  to 
define  the  real  fragment  coordinates  (i.e.  x,  y  and  z).  This  step  has  been  explained  in 
the  previous  chapter  for  station  1 .  Once  this  has  been  done  for  both  stations,  the  issue 
of  matching  fragments  from  station  1  to  station  2  can  be  solved  by  computing  the 

distance  (Dm)  between  the  real  fragment  position  at  station  2  (P2)  and  a  vector  AP^  . 

Vector  AP^  is  defined  as  the  line  between  the  real  fragment  position  in  space  at 

station  1  (PI)  and  the  rear  impact  point  on  target  (A).  Vector  AP2  is  defined  as  the 
line  between  the  real  fragment  position  in  space  at  station  2  (P2)  and  the  rear  impact 
point  on  target  (A).  Figure  7  presents  a  diagram  of  vectors  and  both  x-ray  stations. 
Equations  11,12  and  13  give  the  point  coordinates  while  Equations  14  and  15  are  the 
vectors  joining  those  points. 


Station  1 


Station  2 


A 


Target  center 


Figure  7.  Fragment  position  versus  x-ray  stations 


X, 


(11) 


z 


Where  P 1  is  the  fragment  position  at  x-ray  station  1 . 
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P2  = 


(12) 


X2 

J2 


Where  P2  is  the  fragment  position  at  x-ray  station  2. 


^  =  (-^,0,0)  (13) 

Where  A  is  the  coordinate  for  the  target  center,  b  is  the  distance  between  impact  point 
behind  the  target  and  the  reference  coordinate  system  (x  distance). 

AP^={x^+b,y^,z^)  (14) 

Where  AP^  is  the  line  joining  the  target  center  to  the  fragment  position  at  x-ray 
station  1. 


AP^={x^+b,y^,z^)  (15) 


Where  AP^  is  the  line  joining  the  target  center  to  the  fragment  position  at  x-ray 
station  2. 

The  smallest  distance  ( )  between  the  real  fragment  position  (P2)  and  the  vector 

APy  is  required  as  a  criterion  for  matching  a  fragment  from  station  one  to  a  fragment 
at  station  2.  Obviously,  the  matching  process  involved  looking  at  all  fragments 
combination  to  obtain  the  smallest  distance  Dm  (i.e.  best  match).  The  distance  is 

the  length  of  vector  AP^  minus  PAP2  .  Equation  1 6  presents  how  to  calculate 
PAP2  that  is  the  projection  of  vector  AP2  on  vector  AP^  . 


AP,  (16) 

Equation  1 7  presents  the  distance  (D^)  that  will  be  used  as  a  matching  criterion. 

B. = ^(A]y-p7]yy+(A]%^-PA]yy+(7]y-p7py  (17) 

A  distance  ( )  equal  to  zero  would  mean  that  the  two  fragment  coordinates  (P 1  and 
P2)  are  lying  on  a  line  having  the  origin  at  the  middle  of  target  (A).  Once  the  minimal 


PAP2  = 


AP2»AP, 


AP, 
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distance  is  found,  even  if  it  is  not  equal  to  zero,  the  two  fragments  position  (P 1  and 
P2)  are  identified  as  best  match  and  it  is  therefore  possible  to  obtain  the  distance 
travelled  by  the  fragment  with  real  coordinates.  This  distance  ( )  is  presented  in 
Equation  18. 


D,eal  =  V(^1  -^2)"  +(Ti  -T2)"  +(^i (18) 

Where  Xi,  yi  and  Zi  are  the  real  fragment  coordinates  at  the  X-ray  station  1  (PI)  and 
X2,  yi  and  Z2  are  the  real  fragment  coordinates  at  the  X-ray  station  2  (P2). 

Once  the  real  distance  is  calculated  and  the  time  required  by  the  fragment  to  travel 
from  one  station  to  the  other  recorded,  the  fragment  real  velocity  (V)  can  be  calculated 
with  Equation  19. 


V  = 


D 


real 


t 


(19) 


Where  t  is  the  time  it  takes  the  projectile  to  go  from  one  station  to  another. 
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5.  Impact  point  on  witness  piate 


The  objective  of  this  section  is  to  find  where  the  projectile  impacts  on  the  witness 
plate.  As  seen  in  Figure  8,  the  witness  plate  is  located  according  to  STANAG  4190  [6]. 


A  first  step  is  to  define  a  line  ( Z)3  )  passing  by  the  target  rear  impact  point  and 
reaching  the  fragment  center  of  impact  on  the  witness  plate  (see  Figure  6).  Equation  20 
presents  this  line  equation  according  to  the  reference  coordinate  system.  PI  and  P2  are 
two  points  defined  as  fragment  position  at  station  1  and  2,  respectively,  and  a 
parameter  called  u„,  is  needed  to  obtain  the  interception  on  the  witness  plate.  Flowever, 
before  finding  the  value  of  parameter  u„,  the  witness  plate  normal  must  be  found  with 
Equation  2 1 .  This  is  done  by  performing  the  vectorial  product  of  the  axes  formed  with 
the  point  randomly  chosen  on  the  witness  plate  (i.e.  a,  (3,  y ).  These  three  points  are 
taken  on  the  witness  plate  as  shown  in  Figure  9.  Equations  22,  23  and  24  present  the 
coordinate  of  these  three  points  in  the  reference  coordinate  system.  The  constants 
required  for  these  equations  are  defined  in  Table  1. 


D3  =  P\  +  u„{P2-P\)  = 


X3 

T3 


(20) 
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(21) 


N  = 


-ir.  -«z)(/^.  -«.)  +  (rz  -«.)(A  -«z) 
-  «z  )(r.v  -«>.)-  (A'  -  )(r.  -  «. ) 


P3' 

Z’ 

1 

- 

'  Center 

o  y 

p  • 

> 

a 

Figure  9.  Witness  plate  reference  axes 


a  = 


f  \ 

j 


0 


-Z)AO 


cos 

V  2  y 


f  \ 

j 


0 


cos— 

V  2  y 


Y 


— -  tan  — ,0,1 
0  2 
cos  — 

V  2  y 


(22) 


(23) 


(24) 


Once  the  normal  has  been  defined,  Equation  25  can  be  solved  to  find  the  value  of 
‘  ‘u,„  ’  ’  since  the  scalar  product  of  perpendicular  vectors  is  equal  to  zero.  Equation  26 
presents  the  solution  for  ”  ,  which  can  now  be  used  in  Equation  20  to  find  P3. 


16 


DRDC  Valcartier  TM  2003-123 


{Dl>-a)*N  =  0 


(25) 


^  ^  («,  -Xy)N^+{ay  -yi)Ny  +{a,-z,)N^ 

(X2  -  Xj  )N^  +  ( J2  -  Ji  )Ny  +  (Z2  -  Zj  )N^ 

Where  Um  is  the  parameter  for  finding  a  point  on  the  witness  plate,  x^,y^, Zj  are  the 
real  fragment  coordinates  for  the  x-ray  station  1  (PI),  X2 , 72  ?  ^2  fragment 

coordinates  for  the  x-ray  station  2  (P2)  and  N ^ ,  N ^ ,  N ^  are  the  coordinate  for  the 
vector  orthogonal  to  the  witness  plate. 

Since  the  fragment  position  P3  on  the  witness  plate  is  given  according  to  the  reference 
coordinate  system  ( X ),  a  change  is  required  to  give  the  fragment  position 

according  to  a  coordinate  system  located  on  the  witness  plate  ( x’ ,  ,  z’ )  .  Figure  8 

presents  these  new  reference  axes  with  the  witness  plate  center  as  the  origin  and 
Equation  27  presents  the  fragment  position  on  the  witness  plate  according  to  this  new 
reference. 


py= 


0  .1 

"  r  0  .\ 

r . 

X,  cos - / 

3  2  ■'I 

cos  — 

+  bcos—+  / 

sin  — 

^  2y 

2  1 

73 

^  2j 

.  0 

-  z,  sin  — 
2 


.  e  .  .  e  e  ^  .  e  e  .  e  .  e 

X,  sin - /sin  — cos — hnsin — hz,  cos - /cos  — sin  — 

^2  2  2  2^  2  22 


X3 

73 


(27) 


Where  P3'  is  the  fragment  position  on  the  witness  plate  from  the  center  of  the  plate, 

X3 ,  J3 ,  Z3  give  the  fragment  position  according  to  the  reference  coordinate  system, '  j  ' 

is  the  distance  between  the  witness  pack  center  and  the  impact  point  behind  the  target, 

'  b  '  is  the  distance  between  the  impact  point  behind  the  target  and  the  reference 
coordinate  system  and  finally  '  6^ '  is  the  target  degree  of  inclination. 
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6.  Estimate  of  the  mass 


The  objective  of  this  section  is  to  find  a  way  to  estimate  fragment  mass  based  on  the 
projected  fragment  area  and  the  hole  area  caused  by  the  fragment  on  the  witness  plate. 
Indeed,  there  are  two  horizontal  and  two  verticals  views  representing  the  fragment- 
projected  area.  Therefore,  a  technique  for  measuring  the  projected  area  based  on 
icosahedron’s  method  is  used  according  to  STANAG  4589  [7].  Also  according  to 
STANAG,  the  projected  area  is  the  area  of  a  tumbling  fragment  in  flight  at  a  given 
time  and  rotation.  An  icosahedron  is  a  20-sided  figure,  which  has  10  different  faces. 
The  projected  area  will  then  be  a  mean  value  of  those  corrected  areas.  Since,  only  4 
different  faces  are  available  on  the  film  for  this  present  study,  the  mean  is  taken  on 
those  4  corrected  areas  as  computed  in  Equation  28.  The  fragment  area  found  on  the 
films  will  not  be  used  directly  since  a  fragment  image  is  projected  on  the  film,  which 
is  not  the  real  fragment  area.  A  correction  is  therefore  required  to  obtain  the  real 
projected  fragment  area.  Equation  29  presents  a  way  to  assess  this  correction  for  a 
fragment  seen  at  the  first  station  on  the  horizontal  film  (Figure  1).  A  shape  factor  is 
also  required  to  estimate  fragment  mass.  This  factor  is  presented  in  Equation  30. 


A 


C 


+A„.„  +A„^„  +  A 


^clH 


^clV 


c2H 


^c2V 


4 


(28) 


Where  Ac  is  the  average  corrected  projected  fragment  areas  and^i^i^ ,  ,  A^2h  -> 

A^2y  the  corrected  projected  fragment  area  of  the  4  views  on  the  film. 


A 


clH 


=  A\^ 


f  a-Zi  Y 


(29) 


Where  AIh  is  the  fragment-projected  area  as  seen  on  the  horizontal  film  at  the  first  x- 
ray  station,  Zj  is  one  of  the  coordinates  of  the  fragment  (PI)  and  the  other  parameters 
are  defined  in  Table  1. 


m  = 


/  \3/2 


L 


yjk  j 


Where  m  is  the  estimate  fragment  mass  and  /).  the  shape  factor. 


(30) 
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7.  Procedure 


The  objective  of  this  section  is  to  resume  the  procedure  in  order  to  obtain  the  real 

fragment  velocity  and  its  estimated  mass.  This  is  done  in  step  1  to  5. 

1 .  Select  a  fragment  on  a  view  on  the  first  station  and  then  match  this  fragment  with 
all  fragments  on  the  other  view  and,  for  each  pair  of  fragments,  compute  the 
minimum  distance  between  the  two  vectors  using  Equation  1 .  The  best  result 
should  indicate  a  perfect  match.  If  more  than  one  fragment  give  a  very  good  result, 
two  additional  criteria  could  be  considered.  First,  the  length  of  the  fragment  in  the 
“x”  direction  on  each  view  should  be  similar.  Moreover;  the  distance  from  the 
reference  coordinate  system  to  the  fragment  centroid  should  also  be  very  similar. 
These  last  two  criteria  could  save  computing  time  and  avoid  bad  matching. 

2.  Repeat  step  1  for  station  2. 

3.  The  objective  of  step  3  is  to  associate  fragment  from  a  x-ray  station  to  the  other  by 
computing  the  distance  between  the  fragment  at  station  2  and  the  line  which  joins 
the  fragment  at  station  1  and  the  rear  impact  point  on  target  (Equation  17).  The 
smallest  distance  should  indicate  the  best  match. 

4.  Once  all  fragments  are  matched,  the  results  can  be  validated  using  the  supposed 
impact  point  on  the  WP  (Equation  27).  However,  if  there  is  no  fragment  hole 
located  on  the  WP  at  a  supposed  impact  point,  then  the  match  was  not  right  or  the 
fragment  did  not  have  enough  energy  to  create  a  hole.  Finally,  a  graphic  of  the 
fragment  velocity  calculated  versus  fragment  position  in  space  could  help  to 
confirm  the  matching  accuracy  [1]. 

5.  Once  these  steps  are  completed.  Equation  18  allows  computing  the  real  distance 
travelled  by  the  fragment  between  two  x-ray  stations  and  Equation  1 9  allows 
computing  the  fragment  velocity.  The  last  step  is  to  estimate  fragment  mass  from  4 
fragment  views  on  x-ray  films.  This  is  done  by  first  computing  the  mean  corrected 
fragment  projected  area  with  Equation  28  and  then  applying  Equation  30  to 
estimate  fragment  mass.  The  fragment  area  on  WP  could  also  be  used  has 
additional  information  to  find  the  average  area.  However,  one  should  consider  that 
the  area  seen  on  a  witness  pack  is  bigger  than  the  real  fragment  area  [8]. 
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8.  Conclusion 


This  memorandum  presented  a  way  of  calculating  a  fragment  velocity,  trajectory  and 
estimating  its  mass  from  x-ray  photographs  taken  in  a  physical  installation,  including  a 
target  and  a  witness  pack.  In  order  to  easily  collect  information,  mathematical 
equations  had  to  be  developed.  Using  these  equations  the  association  of  fragment  from 
the  horizontal  film  to  the  vertical  film  and  from  station  one  to  station  two  is  now 
feasible.  Also,  the  results  obtained  with  witness  packs  can  now  be  compared  with 
those  obtained  from  x-ray  systems. 

The  major  objective  of  this  memorandum  was  to  develop  a  method  for  computing  the 
real  fragment  velocity.  An  equation  was  developed  to  allow  this  computation.  Of  equal 
concern  was  to  estimate  fragment  mass  with  the  information  available.  The 
photographs  showed  the  fragment  projections  at  different  times.  Therefore,  two 
overheads  and  two  side  views  were  used  to  find  an  equation  for  the  fragment  mass 
estimation.  However,  there  are  some  limitations  to  the  method  presented  in  this 
memorandum.  Indeed,  since  the  references  are  taken  according  to  the  centroid,  the 
method  is  more  accurate  for  small  projectiles.  Also,  smaller  is  the  projectile  better  is 
estimated  the  fragment  ejection  point.  The  number  of  fragments  must  be  restricted  to 
associate  them  accurately  on  the  x-ray  films.  And  finally,  it  is  assumed  that  fragments 
travelled  along  a  straight  line. 
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9.  Future  work 


Metallic  spheres  have  been  suspended  in  the  impact  chamber  at  different  positions  to 
validate  the  fragment  matching  criteria  for  horizontal  and  vertical  films.  Also,  some 
films  obtained  from  experimental  tests  have  been  analysed  to  prove  the  concept  and 
again  validate  the  process.  Future  work  should  involve  an  experimental  test  where  the 
fragment  velocity  calculated  with  x-ray  would  be  compared  with  the  velocity  obtained 
using  an  independent  system  (e.g.  a  camera  to  record  the  flash  on  a  metallic  witness 
sheet).  Finally,  fragments  could  also  be  recovered  to  measure  the  real  mass. 
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List  of  symbols/abbreviations/acronyms/initialisms 


DND 

BA 

a 

g 

c 

e 


b 


i 

j 

e 

^min 


Dfn 


D 


real 


Ac 

m 

h 

WP 

m 

m 


Department  of  National  Defence 
Behind- Armour 

Vertical  distance  between  x-ray  station  and  fire  line 

Vertical  distance  between  fire  line  and  horizontal  film 

Horizontal  distance  between  fire  line  and  vertical  film 

Horizontal  distance  between  horizontal  x-ray  station  and  fire 
line 

Horizontal  distance  between  impact  point  behind  target  and 
indicator  (x  distance) 

Horizontal  distance  between  two  stations 

Horizontal  distance  between  impact  point  behind  target  and 
WP  center 

Target  degree  of  inclination 

Minimum  distance  between  two  vectors  D\  and  D2 
Minimum  distance  between  a  point  P2  and  the  vector  AP^ 

Real  distance  travelled  by  the  fragment  between  two  x-ray 
stations 

Average  corrected  projected  fragment  area 
Estimated  fragment  mass 
Shape  factor 
Witness  pack 

Parametric  equation  of  a  line  with  extremities  at  Ryi  and  Pin. 
Parametric  equation  of  a  line  with  extremities  at  Rhi  and  Ply. 
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Parametric  equation  of  a  line  going  through  PI  and  P2 


D3 

R 

R 

PIh 


K1 


Plx 


V\,  V2 
A 

PI 


P2 


P3 

AP, 


AP^ 

PAP, 

V 

t 


X 


FH  ’ 


Position  of  the  first  X-ray  station  for  the  vertical  view 

Position  of  the  first  X-ray  station  for  the  horizontal  view 

Position  of  the  projected  fragment  image  on  the  horizontal 
film  at  station  1 

Position  of  the  projected  fragment  image  on  the  vertical  film 
at  station  1 

Scalar  value  to  define  a  specific  position  on  a  line 
Scalar  value  to  define  a  specific  position  on  a  line 
Unitary  vectors 

Target  center 

Real  fragment  coordinate  at  x-ray  station  1  in  the  reference 
coordinate  system 

Real  fragment  coordinate  at  x-ray  station  2  in  the  reference 
coordinate  system 

Real  fragment  coordinate  on  the  WP  in  the  reference 
coordinate  system 

Vector  joining  the  target  center  and  fragment  position 
(station  1) 

Vector  joining  the  target  center  and  fragment  position 
(station  2) 

AP2  vector  projected  on  API  vector 


Fragment  velocity 

Time  taken  by  the  fragment  to  travel  from  station  1  to  station 
2 

Fragment  centroid  position  on  the  horizontal  film  according 
to  X  axe 
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y  FH  ’ 


X 


FV  ’ 


y  Fv  ’ 


a,/3,r 


N 


Fragment  centroid  position  on  the  horizontal  film  according 
to  y  axe 

Fragment  centroid  position  on  the  vertical  film  according  to  x 
axe 

Fragment  centroid  position  on  the  vertical  film  according  to  y 
axe 

Points  on  the  witness  plate 
Orthogonal  vector  to  the  witness  plate 


Scalar  value  to  find  a  point  (/’3  )  on  the  witness  plate  lying 
on  the  line 

Fragment  position  on  the  witness  plate  (coordinate  system 
located  at  the  center  of  the  WP) 
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